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Abstract. We study the joint distribution of descents and inverse descents over the set of 
permutations of n letters. Gessel conjectured that the two-variable generating function of 
this distribution can be expanded in a given basis with nonnegative integer coefficients. We 
investigate the action of the Eulerian operators that give the recurrence for these generating 
functions. As a result we devise a recurrence for the coefficients in question but are unable 
to settle the conjecture. 

We examine generalizations of the conjecture and obtain a type B analog of the recurrence 
satisfied by the two- variable generating function. We also exhibit some connections to cyclic 
descents and cyclic inverse descents. Finally, we propose a combinatorial model for the joint 
distribution in terms of statistics on inversion sequences. 



1. Introduction 



Let & n denote the set of permutations of {1, . . . ,n}. The number of descents in a per- 
mutation 7r = TV i . . . 7r n is defined as des(7r) = \{i : 7Tj > 7r i+1 }|. Our object of study is the 
two-variable generating function of descents and inverse descents: 

A n (s,t)= s dcs ^ +1 t dcs ^ +1 . 
nee n 

The specialization of this polynomial to a single variable reduces to the classical Eulerian 
polynomial: 

n 

A n (t) = A n (i,t)=j2t dcs(7r)+1 = Y,Q tk - 

tt£&„ k=l 

Eulerian polynomials and their coefficients play an important role (not only) in enumer- 
ative co mbinator i cs. T h e classical, univariate poly n omial s are quiet well-studied — see, for 
example, ICarlitzl (Il959h : iFoata and Schutzenbergerl fll970h . This cannot be said for the bi- 
variate generating function for the pair of statistics (des, ides). Here and throughout this 
note we will use the shorthand ides(7r) = des(7r _1 ). Our main motivation to study t hese poly- 
nomia ls is the following conjecture o f Gessel which appeared in a recent article by iBranden 
( 2008 ); see also a nice exposition by Petersen ( 2012[ ). 

Conjecture 1.1 (Gessel). For all n > 1, 

s,t) = 5>M,i(^( s + W + st) n+1 ^- 21 , 



.4, 



where 7ni j are nonnegative integers for all i,j G N. 
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If true, this decomposition would refine the following classical result, the ^-nonnegativity 
for the Eulerian polynomials A n (t). 

Theorem 1.1 (Theoreme 5.6 of Foata and Schiitzenberger ( 1970h ) . 

\n/2] 
i=l 

where are nonnegative integers for all % £ N. 

Before giving their proof, let us recall the recurrence satisfied by the Eulerian polynomials: 

d 

(1) A n (t) = ntA n ^(t) + t(l - t)—A n ^(t) , forn > 2 , 

at 

w ith initial value A-\ (t) = t. 



Foata and Schutzenbergerl (Il970l . Chapitre V) give a purely algebraic proof of Theorem ll.il 



by considering the homogenized Eulerian polynomial, of degree n + 1, 

A n (t-y) = y n+l A n (t/y) 

(2) __ ^dcs(7r)+l^asc(7r)+l 

where asc(7r) denotes the number of ascents (7^ < vr i+1 ) in the permutation n = -K\ . . . 7r n . 
Note that this polynomial is different from and therefore should not be confused with A n (s, t). 
To avoid confusion we use a semicolon and different variables. We include their proof next, 
as we will be applying the same idea to the joint generating polynomial of descents and 
inverse descents in Sectional 

Proof of Theorem The homogenized Eulerian polynomials defined in (jSJ) satisfy the re- 
currence 

(3) A n (t; y) = ty l—A^t; y) + — A n - X (t] y) , for n > 2 , 

which follows from observing the effect on the number of descents and ascents of inserting 
the letter n into a permutation of {1, . . . , n — 1}. Compare this with the recurrence in ([T|). 
It is clear from symmetry observations that A n (t; y) can be written (uniquely) in the basis 

{(m + y) n+1 - 2i } i = 1 .,. w 

with some coefficients 7 nj j. To show that 7^ are in fact nonnegative integers consider the 
action of the operator T = ty (d/dt + d/dy) on a basis element. Apply T on the ith basis 
element we get that 

T[{ty)\t + y) n+1 - 2t ] = i{ty)\t + y) n+2 ^ + 2(n + 1 - 2i){t y y + \t + y) n ~ 2 \ 

which in turn implies the following recurrence on the coefficients: 

(4) 7n+M = hn,i + 2(n + 3 - 2i)7„ )i _i. 

The statement of Theorem 11.11 now follows, since the initial values are nonnegative integers, 
in particular, 7^1 = 1 and 71^ = for i ^ 1. Furthermore, the constraint 1 < i < [|] assures 
that both positivity and integrality are preserved by recurrence (Hj). 

□ 
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Remark. The study of these so-called Eulerian operat ors goes back t o Carlitz as was pointed 
out to the author by I. Gessel. See, for example, Carlitzl ( 1973 ) for a slightly different 
variant of T. Also, the operator t (n + (1 — t)(d/dt)) is closely related to a special case 
of a generalized derivative opera tor already stu died by Laguerre, called emanant or polar 
derivative; see, for example, 6. in iMardenl ( 119351 ). 



Fin ally, we must also mention the "valley-hopping" proof of Theorem ll.ll by lShapiro. Woan. and Getu 
( 119831 . Proposition 4) which is a beautiful construction that proves that the coefficients 7 nj j 
are not only nonnegative integers but that they are, in fact, cardinalities of certain equiva- 
lence classes of permutations. Their proof is part of a more gene ral phenomenon, a n acti on 
of transformation groups on the symmetric group & n studied by iFoata and Strehll ( 119741 ) . 



2. Symmetries of A n (s,t) and a homogeneous recurrence 



The polynomials A n (s,t) were first studied by lCarlitz. Roselle. and Scovillel ( 119661 ). They 
proved a recurrence for the coefficients of A n (s,t) (see equation (7.8) in their article — note 
there is an obvious typo in the last row of the equation, cf. equation (7.7) in the same 
article). The recurrence they provide for the coefficients is equivalent to the following one 
for the generating functions. 

Theorem 2.1 (Equation (9) of lPeterseni (12012^ . For n > 2, 

nA n (s,t) = (n 2 st + (n- 1)(1 - s)(l - £)) K-iM 

d d 
+ nst{l - s)— A,-i(s, t) + nst{\ - t)—A n ^{s, t) 

d 2 

+ st(l-s)(l-t)—A n ^(s,t), 

with initial value Ai(s,t) = st. 

At first glance, this recurrence might not seem very useful at all. However, if we introduce 
additional variables — to count ascents (asc) and inverse ascents (iasc) — we obtain a more 
transparent recurrence. So, let us first define 



(5) 
(6) 
(7) 



A (s t' X y) — ^ ^ s ides(7r)+l^des(7r)+l a ,iasc(7r)+l^asc(7r)+l 



7TS6„ 

E ides(-7r)+l j.des(7r)+l n— ides(7r) ri— des(-7r) 
boxy 

(xy) n+1 A n (s/x,t/y). 



Proposition 2.1. A n (s, t; x, y) is homogeneous of degree In + 2 and is invariant under the 
action of the Klein 4-group V = (id, (12)(34), (13)(24), (14)(23)) ; where the action of a 6 V 
on A n (s,t;x,y) is permutation of the variables accordingly (e.g., o = (13)(24) swaps x with 
s and y with t, simultaneously). 

Proof. The homogeneity is immediate from the second line of the equation above. The 
invariance is a consequence of the symm etry prop e rties o f A n (s, t), such as A n (s, t) = A n (t, s); 
see, for example, equations (12-14) in (Petersen, 20121 ). Note that, due to the introduction 
of the new variables, for n > 4, the polynomial A n (s,t;x,y) is not symmetric. □ 
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Now we are in position to give our homogeneous recurrence. 
Theorem 2.2. For n > 2, 

nA n (s, t; x, y) =(n - l)(s - x)(t - y)A n _ 1 (s, t; x, y) 

(8) / d d\ ( d d\ A , 

+ stxy {dS + d-x) \dt + d-y) A ^ s ^y) 

with initial value Ax(s, t; x, y) = stxy. 

Proof. Consider the bivariate recurrence given in Theorem 12.11 and observe that it can be 
rewritten as 

nA n (s,t) = Un-l){l-s){l-t) + st U + (l-s)JM U + (l-t)-|jJ 

Now we can make both sides of the equation homogeneous using ([7j). Since the two Eulerian 
operators act on different variables each of them can be replaced by their symmetric two- 
variable homogenized counterpart and the theorem follows. □ 

Remark. The invariance of A n (s, t; x, y) under the Klein-group action also follows easily from 
recurrence jS]) directly. Clearly, A±(s, t; x, y) = stxy is invariant under the action of the group 
(in fact, it is symmetric) and the operator acting on A n (s,t;x,y) denoted by 

{ d d \ I d d 

(9) T n = n{s - x)(t - y) + stxy — + 



ds dx J \ dt dy t 
itself is invariant under the action of the Klein-group. 

Finally, Theorem 12 . 2 1 allows us to give a (homogenized) restatement of Gessel's conjecture: 

Conjecture 2.1. 

A n (s, t; x,y)=Y^ ln,ij(stxy)\st + xy) j {tx + sy) n+1 ~ 2i - j , 

where 7 n ,ij G N for all i,j G N. 

For example, we have (cf. page 18 of Petersen ( 20121 )): 

Ax(s,t; x,y) = stxy 

A 2 (s, t] x, y) = stxy(st + xy) 

A 3 (s, t; x, y) = stxy(st + xy) 2 + 2(stxy) 2 

Ai(s, t; x, y) = stxy(st + xy) 3 + 7(stxy) 2 (st + xy) + (stxy) 2 (tx + sy) 

A 5 (s, t; x, y) = stxy(st + xy) 4 + ld(stxy) 2 (st + xy) 2 + 6(stxy) 2 (st + xy)(tx + sy) + 16(stxy) 3 . 

Remark. It is not too hard to see that Theorem 12.21 is. in fact, equivalent to Theorem 12.11 
At the same time, the symmetric nature of the homogeneous operator is more suggestive 
to combinatorial interpretation. It would be nice to find such an interpretation (perhaps in 
terms of non-attacking rook placements on a rectangular board). 
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3. A RECURRENCE FOR THE COEFFICIENTS J n)i j 

Following the ideas of iFoata and Schiitzenberger (I1970L Chapitre V) that were used to 
devise a recurrence for j ni , we apply the operator T n to the basis elements to obtain a 
recurrence for the coefficients J n ,i,j- As a result, we obtain the following recurrence. 

Theorem 3.1. Let n > 1. For all i > 1 and j > 0, we have 

(n + l)7n+i,i,j = (n + i(n + 2 - i - j))n/ n ,i,j-i + + j) ~ n )ln,i,j 
+ (n + A-2i-j)(n + 3-2i- j)Tn,i-ij-i+ 
+ (n + 2z + + 3 - 2% - j)j n ,i-i,j 
+ (j + l)(2n + 2 - j)7n,i-i,j+i + + l)(j + 2 )7n,i-i,i+2 , 

wit/i 71,1,0 = 1; 7i,y = (unless i = 1 and j = 0j and 7 n ,i,j = «/i < 1 or j < 0. 

Proof. Denote the basis elements by B^f = (stxyy(st+xyy(tx+sy) n+1 ~ 2 ' l ~ : > for convenience, 



and recall the definition of T n given in (|9j). 
A quick calculation shows that 

(11) n(s-x)(t- y)B% = n (b%$ - 

To calculate the action of the differential operators on the basis elements, we use the 
product rule, which for second-order partial derivatives is given by the following formula: 

dzw(fgh) = d zw (f)gh + d z (f)d w (g)h + d z (f)gd w (h) 

+ d w (f)dz(9)h + fd zw (g)h + fd z {g)d w {h) 

+ d w (f)9dz(h) + fd w {g)d z {h) + fgd zw (h) , 

where f,g, h are functions, d z = d/dz and d w = d/dw denote the partial differential opera- 
tors with respect to z and w, and d zw = d z d w is the second-order differential operator. 
After some calculations, this gives the following: 

(12) stxy {& t + m) Bi? = ** + 1 - 4 - ^ + ^ + 3 - 



? (n+l) 



+ ( n + 1 _2i-j)(n-2i-j)B\ +hj+1 . 



stx (— —\ 

(13) ' \dsdy dtdx J 



(n+l) 



( n + l_2z-j)(n + 2 + 2z+j)£^ 



j 



Summing (II ip . (1121) and (TTH]) we arrive at the following expression. 
T n [B$] ={n + i(n + l-i- j))B%$ + + j) - n)sg +1) 

+ ( n + i _ 2i - j)(n - 2z - j)^+y +1 + (n + 2 + 2% + j)(n + 1 - 2i - j)B { ^ 
+ j(2n + 3 - j)^-! + J'O " l)^S-2- 

Finally, collecting together all terms T n [S^] which contribute to we obtain ( TlOl) . 

□ 

Remark. If we sum up both sides of (imp for all possible j then we get (jl]) back. 
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One could study the generating function 

G(u,v,w) = ^2ln,i,jU n 



v t w 3 



h3 



with coefficients satisfying the above recurrence. Gessel's conjecture is equivalent to say- 
ing that its coefficients are nonnegative integers. Unfortunately, these properties are not 
immediate from the recurrence ([TO]) . 

4. Generalizations of the conjecture 



Gessell (120121 ) noted that the following equality of iCarlitz et al.l ( 119661 ) 



E 

i,j=0 

can be generalized as follows 



ij + n — 1 
n 



sH j 



A n (s,t) 



[1 - s) n+1 (l -t) n+1 



Let t e& n with des(r) = k - 1. Define A { n\t) by 



E 

i,j=0 



ij + n — k 



n 



Ag\s,t) 



(1 - s)" +1 (l -t) n+1 



Then the coefficient of s l P in is the number o f pairs of permutations (jr, er) such that 
no = r, des(7r) = i and des(a) = j. I Gessell (|2012l ) also pointed out that these polynomi- 
als arise implicitly in iMielnik and Plebahskil ( jl970h : compare (11.10) there with the above 
equation. 

This suggests that Conjecture 11.11 h olds in a more g eneral form (this version of the con- 
jecture appeared as Conjecture 10.2 in iBrandenl (120081 )). 

Conjecture 4.1 (Gessel). Let r e & n . Then 

s dcs w +1 t des ^ lr ) +1 = J2in,i/sty( s + ty(i + st) n+1 -i- 21 , 



hJ 



where j^ij are nonnegative integers for all i,j G N. Furthermore, the coefficients do 
not depend on the actual permutation r, only on the number of descents in r. 

In the special case when r = n(n — 1) ... 21 (and hence des(r) = n — 1) the roles of descents 
and ascents interchange. 

Theorem 4.1. For n > 2, 

nA^(s,t;x,y) ={n-l){x- s){t - y)A^\s,t;x,y) 



(14) 



with initial value A^\s,t; x,y) = stxy. 

In particular, we have the following identity. 
Corollary. 



A%\s,t;x,y) = A n (s,y; x,t). 
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4.1. A type B analog. iGessell ( 120121 ) also noted that there is an analogous definition for the 
hyperoctahedral group Q3 n . The elements of 23 n can be thought of as signed permutations of 
{1, . . . , n}, and the type B descents are defined as des# (cr) = {i G {0, 1, . . . , n} : o~i > cr i+1 } 
with cr := for o = o\ . . . o n G *B„. 



E 

i,j=0 



2ij + i + j + l + n — k 



n 



s¥ 



#(s,t) 



[1 - s)"+l(l - t) n+1 



where 



B^(s,t) = s dcSBia H deSB{a ' lT \ 

cre<B n 

with r G 23 n such that des^r) = k — 1 (here des# denote s the des c ents o f type B). 

Therefore, mimicking the proof of Theorem 12.11 given bv lPetersenl (l2012h . we get an analog 

of Theorem 12.11 for the type B two-sided Eulerian polynomials, B n (s,t) = Bn (s,t) 
Theorem 4.2. For n > 2, 

nB n (s,t) = {2n 2 st — nst + n)B n _i(s,t) 

d 

+ (2nst(l -s) + s(l - s (1 - ty—Bn^t) 

OS 

(is) a 

+ (2nst(l - t) + t(l - s)(l - t))—B n ^(s,t) 

+ 2st(l-s)(l-t)— Sn^S.t). 

initial value Bi(s,t) = 1 + st. 

Proof. Following the proof for the case of the symmetric group in (IPetersenl . [20121 . eq. (9)), 
we use the corresponding identity of binomial coefficients: 

[2ij + % + j + n\ . (2%j + % + j + n — 1\ (2ij + i + j + n — l 

Multiplying both sides by the monomial s l P and summing over all integers i,j we get 
'2ij + i + j + n N 



E 



n 



2ij + i+ j + n — 1 
n — 1 



2zj + + — 1 

?2—l 



i,j=0 



from which we obtain the following recurrence for F n (s,t) = B n (s,t)/(1 — s) n+1 (l — t) n+1 : 



d 2 d d 

nF n (s, t) = 2st—-F n _ 1 (s, t) + s— F n _x(s, t) + t— F n _i(s, t) + nf„_i(s, t) . 
as at as at 
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Now substitute back the expression for F n (s, t), multiply both sides with (1 — s) n+1 (l — t) n+1 
and with a little work we get that 

nB n (s,t) = (2n 2 st + nt(l - s) + ns(l - t) + n(l - s)(l - t)).B n _i(s, t) 

d 

+ (2nst(l -s) + s(l - s)(l - t))—B n ^(s,t) 

d 

+ (2nst(l -t)+ t(l - s)(l - t))— S n _!(s,t) 
+ 2st{l-s){l-t)—B n _ 1 {s,t). 

□ 

It would be of interest to find a homogeneous version of this theorem (an analogue of 
Theorem 12.21) and a recurrence for the corresponding 7 n ,i,j coefficients in the case of type B. 

4.2. Cyclic descents. One can also consider two-sided Eulerian-like polynomials using 
cyclic descents. A cyclic descent of a permutation 7r in & n is defined as 

cdes(7r) = \{i : m > 7T(» + i) mo dn}| = des(7r) + x{^n > 7Ti) , 

where 

1, if a > b, and 
0, otherwise. 

The following theorem refines a result of iFulman fl200l Corollary 1). 
Theorem 4.3. For n > 1, 

(n + 1)A„(S, t) = ^ gCdesCTr-^^deSTr _ 
7re6„ + i 

Lemma 4.1. Let a = 23 ... nl denote the cyclic rotation in <5 n (for n > 2). T/ien 

(cdes(7r), cdes(7r -1 )) = (cdes(7r<7), cdes((7rcr) -1 )). 

In oi/ier words, the cyclic rotation simultaneously preserves the cyclic descent and the cyclic 
inverse descent stastics. 



X{a > b) 



Remark. Lemma [4.11 is essentially the same as Theorem 6.5 in (ILam and Postnikovl . 120121 ). 
We give an elementary proof of it, for the sake of completeness. 

Proof. The part that cdes(7r) = cdes(7r a) is obvious since cyclical rotation does not effect the 
cyclic descent set. For the other part, it is equivalent to show that cdes(7r) = cdes(a _1 7r). In 
other words, the cyclic descent statistic is invariant under the operation when we cyclically 
shift the values of a permutation, i.e., add 1 to each entry modulo n. For ir = tti . . . n n 
an arbitrary permutation in & n denote the entry preceding n and following n by a and b, 
respectively. Then ix = i\\ . . . anb . . . ir n and a _1 7r = (m + 1) . . . (a + 1)1(6 + 1) . . . (vr„ + 1). 
Clearly, in all but one position the cyclic descents are preserved, same is true for the cyclic 
ascents. The a /* n cyclic ascent is replaced by the (a + 1) \ 1 cyclic descent and similarly, 
n \ b gets replaced by 1 /* (b + 1). Thus, the total number of cyclic descents remains the 
same. □ 



Proof of Theorem \4.3[ Using Lemma H~T1 we can apply the cyclic rotation to any permutation 
in & n +i until 7r n+ i = n + 1. This will map exactly n + 1 permutations in (5 n +i to the 
same permutation m . . . vr n (n + 1). Clearly, cdes(7Ti . . . n n (n + 1)) = des(7Ti . . . 7i n ) + 1 and 
cdes((7T! . . . 7r„(n + 1)) _1 ) = des((7Ti . . . 7r n ) _1 ) + 1 and the theorem follows. □ 

5. Connection to inversion sequences 

We conclude by proposing a combinatorial model for the joint distribution of descents and 
inverse descents. 

A permutation 7r e & n can be encoded as its inversion sequence e = (ei, . . . , e n ), where 

e - = \{i : % < j )7ri > 7Tj } | . 
Let I n = {( ej , . . . , e„) £ Z n : < < i — 1} denote the set of inversion sequences for & n . 



Recently. ISavage and Schuster! ( 20121 ) studied the ascent statistic ascj(e) = \{i : < e i+ i}\ 



for inversion sequences (and their generalizations) and showed that this statistic is Eulerian, 
i.e., it is equidistributed with the descent statistic over permutations. We use the subscript / 
to emphasize that this is a statistic for inversion sequences which is different from the ascent 
st atistic for permutations used earlie r in the paper. 

Mantaci and Rakotondrajao ( 200ll ) also studied this representation of permutations under 



the name "subexceedant functions" . They considered the statistic that counts that distinct 
entries in e G I n , dst(V) = |{e$ : 1 < i < n}\. They gave multiple proofs of the following 
observation (which they attributed to Dumont) that this statistic is also Eulerian. 



Proposition 5.1 (Dumont). 



x dst(e) 



In fact, the joint distribution (ascj, dst) over inversion sequences seems to agree with the 
joint distribution (des, ides) of descents and inverse descents over permutations. 

Conjecture 5.1. 

A n (s,t) = ^2 s dst ( e) t asc/(e)+1 . 

e€ln 

This observation clearly deserves a bijective proof. Such a proof might shed light on a 
combinatorial proof of recurrence ©. Note that it is not even clear to begin with why the 
right-hand side should be a symmetric polynomial in variables s and t. 
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